Mesoscopic fluctuations of the ground state spin of a small metal particle 
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We study the statistical distribution of the ground state spin for an ensemble of small metallic 
grains, using a random-matrix toy model. Using the Hartree Fock approximation, we find that 
already for interaction strengths well below the Stoner criterion there is an appreciable probability 
that the ground state has a finite, nonzero spin. Possible relations to experiments are discussed. 

PACS numbers 73.23.-b, 71.10.-w, 71.24+q, 75.10.Lp 



According to Hund's rule,0 electrons in a partially filled 
shell in an atom form a many-body ground state with 
maximum possible spin. The maximum spin is preferred 
because it allows a maximally antisymmetric coordinate 
wavefunction in order to minimize the electrostatic re=j 
pulsion between the electrons. In recent experiments,EI 
Hund's rule was also observed in a cylindrically-shaped 
semiconductor quantum dot, or "artificial atom". The 
close similarity with real atoms is due to the degeneracy 
of single-particle levels, caused by the the high degree of 
symmetry of the device. 

In generic ultrasmall systems such as. -small metal 
grains,oo [Semiconductor quantum dots|m or carbon 
nanotubesLrB there is no systematic degeneracy due to a 
spherical (or cylindrical) symmetric potential. However, 
even in the absence of degeneracies, a nonzero value of 
the ground state spin may occur, as long as the gain 
in electrostatic energy is larger than the loss in kinetic 
energy when an antisymmetric coordinate ground state 
wavefunction is formed. Such a ground state is most 
likely to be detected in ultrasmall metal and semiconduc- 
tor devices, since in those systems, unlike in macroscopic 
samples, the spacing between single-particle energy levels 
and the typical interaction energies can be larger than the 
temperature. In fact, the possibility of such a "weakly 
ferromagnetic" ground state has been suggested as an ex- 
planation for some recent experiments, thatxguld not be 
explained by simple noninteracting models. ana In addi- 
tion, a nonzero ground state spin from numerical simu- 
lations of a few particles in a chaotic, dot ,13 and a theory 
of spin polarization in larger dotscil were already men- 
tioned in the literature. The stability of the zero spin 
ground state in a quantum dot was analyzed for weak 
interactions in Ref. [j| 

In this paper, we consider small metal grains in the 
mesoscopic regime, in which fluctuations of wavefunc- 
tions and energy levels, caused by, e. g., disorder or an 
irregular shape, control the behavior of kinetic and inter- 
action energies at the vicinity of the Fermi energy. As a 
result, the ground state spin becomes subject to sample- 
to-sample fluctuations. Then, the relevant quantity to 
consider is the statistical distribution of the ground state 
spin for an ensemble of small metal grains or chaotic 
quantum dots, rather than the spin of a specific sample. 

Our starting point is a simple toy model that captures 
the essential mechanisms for mesoscopic fluctuations of 



the ground state spin. In second-quantized form, our 
model Hamiltonian TL reads 
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where c£ a (c n(T ) is the creation (annihilation) operator 
for an electron with spin er at site n. The indices m, n are 
summed over M sites. The Hamiltonian TLq contains the 
kinetic energy and the impurity potential. We describe 
the electron-electron interaction by an on-site (Hubbard) 
interaction, uM . While the long-range Coulomb inter^. 
action can be trivially included via a charging energy,cl 
the model ([!]) does not include the Coulomb interac- 
tion at intermediate distances, which leads to arGjaussian 
level spacing distribution at the Fermi energy.113 In this 
work, we report a calculation of the ground state spin 
of the Hamiltonian ([!]) using a restricted version of the 
Hartree- Fock (HF) approximation with a random-matrix 
assumption for the eigenvalues and eigenvectors of the 
self-consistent HF Hamiltonian. 

We first present our main result. It consists of an 
equation that relates the candidate ground state energies 
Eq(s) for different values of the total spin s in terms of 
eigenvalues of a hermitian random matrix with level 
spacing A, the interaction parameter A = it/A, and a 
(nonuniversal) numerical constant c that describes the 
density response to a local perturbation of the impurity 
potential in 7io, 



E G (s)~E G (s ) = J2(el 



N+ti+2s 



r 



AA 



2(a -3 ) 
0(1 - A 2 c 2 ) 



(2) 



The total number of electrons is 2(N + sq), so being or 
1/2. The spin of the true ground state is found by min- 
imizing Eq. (j^) with respect to s. The parameter 0—1 
(2) if time-reversal symmetry is present (absent). The 
effect of spin orbit coupling and Zeeman splitting is not 
included here. (The case = 2 is only relevant for semi- 
conductor quantum dots in a weak magnetic field, that 
affects orbital motion, but causes no Zeeman splitting. 
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FIG. 1. The probability distribution P(s) of the ground 
state spin of a small metal grain, computed from Eq. (^) for 
three different values of the interaction parameter A. The 
upper (lower) histograms are for the presence (absence) of 
time-reversal symmetry. Solid histograms are for integer spin, 
dotted ones for half-integer spin. (The density-response pa- 
rameter c has been set to zero; finite c results in an even 
higher probability to find nonzero spin.) 

It is not relevant for small metal grains,! as laboratory 
magnetic fields do not affect orbital motion in this case.) 
Equation (||) reflects the competition between kinetic en- 
ergy (first term on the r.h.s.), which favors small s and 
the on-site interaction (second term), which favors finite 
s. The interaction term, in turn, consists of two parts: 
A term quadratic in s, which describes the exchange in- 
teraction, and a term linear in s, which describes the 
additional "dressed" Coulomb repulsion of two particles 
with the same spatial wavefunction. For large s, the con- 
tribution from the kinetic term is approximately s 2 A, so 
that for u > A a finite fraction of the total number of 
spins will align, rather than a small number of spins as in 
the case u < A. The instability at u = A is known as the 
Stoner instability. In Ref. |[ a result similar to Eq. (^) 
was obtained for s = 1, but with a different and fluctuat- 
ing interaction term. The difference is due to the absence 
of a self-consistent approximation scheme in Ref. |9[ 

As a consequence of the additional dressed Coulomb 
repulsion of particles with the same wavefunction, we 
find that already for interaction strengths considerably 
below the Stoner instability u = A, there is an appre- 
ciable probability of nonzero ground-state spin. This is 
illustrated in Fig. [I], where the distribution of the ground 
state spin at three different values of the interaction pa- 
rameter A is shown: Already at the quite modest inter- 
action strength u w 0.4A a ground state spin s = 1 is 
more likely than s = 0. 

The effect of a weak magnetic r%ld is twofold: First, 
it changes the statistics of the £° Jl3 and, second, it sup- 
presses the interference in the "Cooper channel" , leading 
to a factor of two reduction of the interaction term linear 
in s [last term in Eq. (Q)]. Both effects favor lower s than 
without a magnetic field. However, even in the absence pf 
a magnetic field we expect that, similar to 3D metals,til 
inclusion of the electron-electron interaction beyond the 
HF approximation will also lead to a suppression of the 



interference in the Cooper channel (logarithmically in the 
system size M), and hence to a prefactor in that term 
that is smaller than 2. 

Let us now turn to the details of our calculation. To 
find the ground state of the Hamiltonian ([!]) we use a 
simplified version of the HF approximation: We assume 
that the ground state has the form of a Slater determi- 
nant of single-particle wave functions Vvi an d VVd °f 
particles which have either spin up or spin down. In this 
case the self-consistent HF equations read 



HnF,a(n, m) = Ho(n, m) + uMp- a (n)5(n, to), 

At 00 =^2U,a\lp,M{n)\ 2 . 
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The occupation number / MjCr is 1 (0) if the level p, a is 
occupied (unoccupied) and S(n, to) is the Kronecker delta 
function. The ground state energy Eq is given by 
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Our strategy is as follows: We start from a reference 
state with zero spin, in which N particles of each spin 
are placed in the same levels = e^.i and with the 
same wavefunctions ip^ — ip^i- We assume, that for 
this symmetric case the eigenvectors and the eigenvalues 
of Hhf arc distributed like those of a random matrix, ex- 
cept that the energy levels below Ep are shifted upwards, 
by a small constant amount, relative to the levels above 
Ep , see Eq. ( |l4| ) below. (If we would have included long- 
range Coulomb interactions via a charging energy, the 
shift would have been much larger and in the opposite 
direction. Omission of the charging energy has no con- 
sequence in our case, as we compare ground states with 
the same number of particles.) The assumption that the 
single-particle eigenvalues and wavefunctions in a self- 
consistent potential for a mesoscopic system below and 
above Ep obey random matrix statistics, even though 
they may be quite different from their counterparts in the 
noninteracting system, was checked numerically fopshort 
range interaction models somewhat similar to oursJ13 The 
energy shifts in our case result from the spin-degeneracy, 
which was not present in these calculations. 

Starting from this reference state, we build other states 
by the subsequent addition and removal of electrons. 
We first discuss the addition of a single up spin in the 
(N + l)st level. The first question that needs to be an- 
swered is how this addition affects the self-consistent den- 
sity p a {n). The density change Sp^(n) consists of a direct 
and an induced contribution, while 5p[(n) has an induced 
density shift only, 

5pi = Sp tdil + 5p Uad , Sp^ dil (n) = I^w+i.tOOI 2 ! 

$Pl = SPl.ind- (5) 

Since the density shifts change the HF Hamiltonians by 
an amount STi.nF,a-(n, m) = uM6p- a (n)5(n,m), we ob- 
tain the following self-consistency equations for 5p a , 
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Both 5pi and 8pi are of order 1/M . In Eq. (^) we have 
computed the induced density change to first order in 
STLyw- Higher order terms do not contribute to <5pi n d,er 
to order 1/M and are neglected. To evaluate Eq. (g), we 
first sum the r.h.s. over the space index m and then over 
the energy levels £ M and e u . Because the eigenfunction 
has a random sign, a single term in the latter summation 
is of order AM~ 3 / 2 A/(e M — e„), which is not relevant in 
the limit M ^> 1, even if and e v are both close to 
the Fermi level Ep. For the summation over all levels we 
may perform an average over the wave functions (since 
the denominator is a slowly varying function of p and v 
away from the Fermi level). This average is done usiag 
that for general p ^ v and in the limit M ^s> 1 one has,E2l 

(^, ff (»)i, ff ( m )C ! ,( m )i,a( 11 )) = ^^p^ ~ ^p- 

Putting everything together, we find the following so- 
lution of the self-consistency equations (|J) and (^), 
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where c is a numerical constant of order unity defined by 
c= hm 2 A[ EF d£l r°° d£2 E(E2kM. (8) 



M^oo M 



£2 - £l 



In this equation, p(e) is the mean density of HF energy 
levels. [The mean level spacing A is taken at the Fermi 
energy, A = l/p(ep).] The constant c gives the linear 
density response Sp(n) to a shift of the impurity poten- 
tial H.o(n, n) at that same site, Sp(n) = c5TLo{n, n)/MA, 
which can be verified using first order perturbation the- 
ory with respect to SHo- Notice that c is not a universal 
constant, but depends on an integration of the density of 
states over the entire bandwidth. For example, for the 
Wigner semicircular density of states we find c = 4/3 if 
the Fermi energy Ep is at the band center and c — > if 
Ep is at a band edge. Equation (jjj) expresses that the in- 
teraction enhances the fluctuations of the spin density: if 
l^jv-fi f(7i)| is larger than average, the on-site repulsion 
reduces pi(n), which in turn causes an increase of Pi(ri), 
and so on. At the same time, the interaction reduces 
fluctuations of the charge density p\ + Pi- 

At Ac = 1, which may occur before the Stoner instabil- 
ity A = 1 if c > 1, the density changes diverge. Although 
this instability signals a breakdown of our approach, it is 
not clear whether it will also lead to a true macroscopic 
ground state spin. Below, we restrict our discussion to 
the case Ac < 1. 



Next we address the HF energy levels e^ a and find 
. , . 2 cA 
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The shift of £jv+i,ct is extra large, since for that level 
the interaction effects are enhanced by the spatial fluc- 
tuations of the wavefunction. Equation (JTJ) is the re- 
sult of first order perturbation theory in Sp a ; second or- 
der perturbation theory 8e^ a gives a correction of order 
A 2 AlnM/M, which we may neglect in the limit M 3> 1. 
In the same way, one finds that the changes in each indi- 
vidual wave functions is not significant for M 3> 1. 

Finally, we consider the change in the ground state 
energy Eq. Because the summation over p in Eq. (Q) 
extends over 0(M) levels, it is important to follow the 
shifts in the HF levels to second order perturbation the- 
ory, although this level of accuracy was not needed for 
the shift of each level individually, cf. Eq. (^|). Putting 
everything together, we find that 
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c 2 A 2 )] -1 . 



(10) 



(No terms proportional to log AT appear here since they 
cancel in the summation over the energy levels.) With 
the help of Eq. (||), we can rewrite Eq. (|Io| ) as 8Eq = 
£jV+i,t + \ 8en+i,], which is the average of the energy for 
the newly occupied level £jv+i,t before and after its oc- 
cupation. This may bejaterpreted as a simple extension 
of Koopmans' theoren£3 to the present case, where the 
modification in each one-electron wavefunction is small 
(of relative order M^ 1 / 2 ), but the resulting contribution 
to 8Eq cannot be neglected to the order we are interested 
in. In the usual form of Koopmans' theorem, where one 
ignores any change in the one-particle wavefunctions, the 
HF energies of the lowest unoccupied state before addi- 
tion of the electron, and of the highest occupied state 
after addition are identical. (The usual Koopmans' the- 
orem is correct for an infinite system, in general, or for a 
translationally invariant finite system, as the one-electron 
states are trivially plane waves in that case.) The sim- 
ple extension of Koopmans' theorem also works, in our 
model, for the addition of several electrons. 

We have repeated these calculations for the addition of 
two electrons with opposite spin in the (N + l)th level, 
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(12) 
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As in the case of the addition of a single particle, the 
individual wavefunctions do not change to order M -1 / 2 . 
Equation ( |l2| ) allows us to find the statistics of the HF en- 
ergy levels s^a- in our reference system with N electrons 
of each spin: The only distribution that is consistent both 
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with the assumption that the e^ j(T obey random-matrix 
statistics away from the Fermi level and with the shifts 
of Eq. (|5I 



is one where the have the form 
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where the e° have random matrix statistics and = 1 
(0) if the level fi is (un)occupied. In other words, the 
distribution of the is the same as that of the eigenval- 
ues of a random matrix, with all occupied levels shifted 
upwards by an amount (2//3)AA/(l + cA). 

With the knowledge we have gained above, there is lit- 
tle work left for the calculation of our main result, Eq. 
(|J). Some remarks about the validity of this result are 
appropriate. First, to make a connection between our 
random matrix toy model and a laboratory made quan- 
tum dot we must identify M = (L/Xp) 3 asrwe expect 
the length scale for wavefunction correlations!!!] and the 
range of the screened Coulomb interaction to be of order 
of the Fermi wavelength A^-t-J Second, while our solution 
is complete within the HF approximation, one must bear 
in mind that this approximation scheme does not include 
correlation effects, such as the Cooper channel renormal- 
ization, as we discussed below Eq. (||) . Those correlation 
effects are not expected to affect our result to first order 
in A, but it can not be excluded that they are important 
in the higher order terms in Eq. (^), which involve the 
factor cA. 

We close this paper with a discussion of the physical 
consequences of a ground state with spin s > 1/2 and of 
the experimental situations in which it can be observed. 

First, the temperature T needs to be smaller than the 
separation of the ground states for different spin s. This 
separation, which is typically smaller than the single par- 
ticle level spacing A, is a fluctuating quantity. Very small 
values are possible, because, unlike in noninteracting ran- 
dom systems, there is no level repulsion if states of dif- 
ferent spin are involved. 

For sufficiently low T, the magnetization of the grain is 
proportional to the spin s of the ground state. However, 
s will also affect other properties which are more easily 
accessible in an experiment, like current-voltage char- 
acteristics: A nonzero ground state spin can serve to ex- 
plain the absence of an even-odd structure in the a H Hit-inn 
spectra of Coulomb-blockaded quantum dots,tr E3Bj Bj or 
the presence of kinks in the parametric dependence of 
Coulomb blockade peak positions, as was noted in Ref. 
E2|. Spin is also relevant for conductance measurements 
at a finite bias voltage, which allow fc>f-# "spectroscopy" 
of the quantum dot or metal grainJj'ta In the presence 
of a magnetic field, the ground state is split by the Zee- 
man energy, and the differential conductance will show 
two peaks, whose relative intensity differs by a factor 
2sjv + i + 1 or 2sjy + 1, depending on whether thej-tjmncl- 
ing onto or from the grain is the faster process.lla Even 
without an external magnetic field the ground state may 
be split, e.g., by spin-orbit coupling or magnetic impu- 
rities, and thus give rise to a multiplet of peaks in the 



differential conductance.^ The peak separation within 
a multiplet is controlled by the strength of the splitting 
mechanism and may be much smaller than A. 
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